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E  LONGA  TIONAL  VISCOSITY  RHEOLOGICAL 


PARAMETER  STUDY 

1.  INTRODUCTION 

The  following  objectives  were  set  under  this  research  contract: 

a.  To  establish  a  hypothesis  relating  the  deformation  of 
Newtonian  and  viscoelastic  liquids  to  an  elongational  viscosity 
parameter 

b.  To  design  and  develop  a  technique  to  determine  this 
parameter  for  a  series  of  liquids  over  an  extended  range  of 
elongation  rates. 

These  objectives  were  aimed  at  elucidating  the  process  of  dis¬ 
persion  of  a  liquid  into  droplets,  and  the  role  that  liquid  properties  play 
in  affecting  that  process. 

For  a  variety  of  circumstances  it  may  be  desirable  to  take  a  coher¬ 
ent  mass  of  liquid  and  cause  it  to  be  distributed  into  a  spray  or  cloud  of  fine 
particles.  The  ultimate  disposition  of  this  cloud,  its  spatial  dispersion, 
settling  velocity,  and  ultimately  its  interaction  with  surrounding  surfaces, 
depends  strongly  upon  the  particle  size  distribution  created  by  the  breakup 
mechanism.  The  particle  size  distribution  depends  upon  the  mechanism 
used  to  form  the  initial  cloud,  and  upon  the  physical  properties  of  the 
liquid  mass.  It  is  well  recognized  that  surface  tension  is  a  primary 
liquid  property  affecting  drop  breakup.  In  viscous  liquids  one  observes 
that  the  liquid  viscosity  exerts  an  effect  on  drop  size.  In  liquids  containing 
polymeric  modifiers  it  has  become  clear  that  an  additional  property,  re¬ 
lated  in  some  way  to  viscoelastic  behavior,  is  necessary  in  order  to 
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describe  braakup  characteristics.  It  is  by  no  means  clear ,  however,  what 
property  one  needs  to  measure.  The  primary  goals  of  this  program  are 
to  clarify  this  last  point. 

A  significant  literature  exists  in  the  general  area  of  drop  breakup 
by  high  velocity  airstreams.  *  ^  All  of  these  studies  suggest  that  drop 
breakup  involves  two  mechanisms:  ligament  formation  associated  with 
pressure  and  velocity  fluctuations,  and  drop  formation  by  Rayleigh  in¬ 
stability  of  the  ligaments. 

Ligament  formation  would  appear  to  be  one  mechanism  that  de¬ 
termines  drop  size.  Thus,  it  is  necessary  to  develop  an  understanding  of 
the  role  played  by  fluid  properties  in  affecting  the  ease  with  which  some 
energetic  mechanisrr  can  "pull"  ligaments  from  the  main  body  of  the  liquid 
mass. 

One  may  regard  ligament  formation  as  a  stretching  flow,  similar 
to  that  ir.  which  a  cylindrical  specimen  of  material  is  elongated  as  a  cylinder 
of  increasing  length  and  decreasing  diameter.  Because  this  is  an  elongational 
flow,  and  not  a  shearing  flow,  it  is  the  elongational  rheology  of  the  liquid 
that  plays  the  major  role  in  determining  the  resistance  of  an  element  of 
liquid  to  ligament  formation. 

In  the  past,  elongational  rheology  studies  have  been  restricted  to 
molten  polymers  or  rubbery  elastomers  of  sufficiently  high  viscosity  (greater 
than  1 0^  poise)  that  can  be  formed  into  a  cylindrical  specimen  and  pulled 
in  a  stress- strain  type  of  test.  The  liquids  of  interest  here  are  of  relatively 
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low  viscosity  (typically  10  to  100  centipoise)  and  cannot  be  forrred  into 
a  specimen.  The  challenging  task  at  hand,  then,  is  the  development  of 
a  test  procedure  which  provides  a  measure  of  the  elongational  rheology 
of  such  liquids. 

We  examine,  first,  some  simple  models  of  drop  deformation 
under  conditions  that  would  suggest  the  importance  of  elongational  viscosity. 
Then  we  describe  the  principles  of  the  technique  proposed  for  the  measure¬ 
ment  of  elongational  response  of  a  liquid.  Progress  toward  creation  of  an 
apparatus  for  performing  this  measurement  is  then  described. 


2.  PREVIOUS  RELEVANT  RESEARCH 

The  most  relevant  studies  are  cited  in  the  list  of  references.  ^ 

The  generally  agreed  upon  picture  that  emerges  is  the  following:  a  mass 

of  liquid  is  first  pulled  into  ligaments  which  break  down  to  large  drops. 
Subsequently,  a  droplet  undergoes  two  important  events.  The  first  is 
acceleration  of  the  drop,  which  may  lead  to  Taylor  instability  of  the  surface. 
It  would  appear  that  a  criterion  of  breakup  through  this  mechanism  depends 
upon  the  bond  number 


L  . 

Bo  =  g  o  r  / (7 
o 

where  g  is  the  drop  acceleration 

p  is  the  drop  density 

r  is  the  drop  radius 
o 

O  is  the  inlerfacial  tension 


(1) 


o 


We  note  that  this  criterion  does  not  involve  any  rheological  parameter. 

The  second  event  of  importance  is  the  deformation  of  the  drop 
into  a  thin  disk.  Studies  cited  suggest  that  the  disk  may  deform  into  a 
"bag"  which  then  disintegrates  into  droplets.  Alternatively,  some  studies 
indicate  that  liquid  is  stripped  from  the  rim  of  the  disk.  In  either  case 
disk  formation  is  central  to  this  mode  of  breakup. 

In  what  follows  we  model  the  process  of  drop  deformation  into  a 
thin  disk,  and  suggest  that  the  presence  of  polymeric  additives  retards 
this  deformation.  We  show  that  disk  formation  is  an  elongational  flow. 

Thus  we  hypothesize  that  elongational  rheology  is  central  to  the  deforma¬ 
tion  process,  and  that  it  is  reasonable  to  seek  a  connection  between  elonga¬ 
tional  rheology  and  drop  breakup. 

3.  A  MODEL  FOR  DISK  DEFORMATION 

To  provide  a  tractable  problem  we  take  the  following  model.  The 
spherical  drop  is  replaced  by  a  cylinder  of  equal  length  and  diameter  (see 
Figure  1).  At  some  initiaL  instant  t  =  0  a  uniform  pressure  is  imposed  on 
the  circular  faces  of  the  disk,  causing  it  to  deform  into  a  thinner  disk. 

We  neglect  the  viscosity  of  the  surrounding  air. 

Newtonian  fluid: 

1  d  ,  z 

Continuity  equation  —  — — (r  v  )  +  — —  =  0  (2) 

r  or  r  oz 

Assume  a  solution  with  properties 
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(3) 


v  r  =  r  f(z) 
p  =  constant  =  P 


(4) 

(5) 


Neglecting  disk  inertia  we  find  the  radial  momentum  equation  that  must  be 
satiaified  to  take  the  form 


=  0 


(6) 


A  solution  satisfying  these  equations  is 


v  =  a  r  (7) 

r 

where  we  have  assumed  that  no  shear  stress  is  exerted  by  the  surrounding 
inviscid  gas: 


3  V 

T  =  0  =  ~ 
rz  3  z 

The  equation  of  continuity  is  integrated  to  yield 


(8) 


2  77  r 


dz 


2 

i  77  r 


H  =  0 


(9) 


from  which  it  follows  that 


a 


-I  IrJi) 

2  H  ' 


(10) 


Here,  H  is  the  half-thickness  of  the  disk  at  any  time  t,  and  (-H)  is  half 
the  rate  of  thinning. 


)2 


1  3 


To  find  C  we  use  the  condition  that  r'  =  r  at  t'  =  t ,  where  H*  =  H; 


this  gives 


ii  _  r  H(t) 

r  jH'(t')  _ 


In  a  similar  manner  we  find 


z_  _  H’(t’) 
z  ‘  H(t) 


The  components  of 


are  found  using 


The  pressure  causing  the  deformation,  P,  may  be  written  as 


P  =  -  T 

z  z 


r 

r  r 


) 


where  we  assume  T  =  0  . 

r  r 

Now  we  must  introduce  some  realistic  viscoelastic  constitutive  equation. 
Related  studies  in  our  laboratories^  suggest  that  a  useful  model  is  the 
integral  model  due  to  Wagner 


(22) 


P  =  r  -  t  =  /  M(t  -  t')  (y;  ,  -  y'  )  dt' 

rr  zz  l  [ojrr  [ojzz J 


(23) 


dt' 


The  simplest  ease  to  examine,  computationally,  and  one  that  suffices 
to  illustrate  the  role  of  elongational  rheology  in  disk  deformation,  seeks  the 
pressure  required  to  achieve  a  specified  rate  of  thinning.  Thus  we  assume 


-H/H  =  €  =  constant 

o 


(24) 


and  find,  upon  integrating  Eq.  23,  that 


V  (T;  A  )  = 


exp 


1  -  2aJ  T 


1  -  2A  / 

-rri  1 1  - exp  ["  (‘ +  A) T] 

We  have  defined  the  following  dimensionless  groups: 

V  -  P(t ;/ n  ? 

o 

A  =  \  i 

o 

J  =  t/A 


(25) 
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;  'i 


Cl 


We  have  chosen  M(t  -  t')  to  take  the  form 

M  =  — exp  [  -  (t  -  t')/X]  (26) 

A^ 

which  corresponds  to  the  codeformational  Maxwell  fluid. ^ 

Figure  2  shows  il  (T,  A)  versus  T.  For  the  Newtonian  fluid, 
obtained  from  Eq.  26  as  the  limit  when  A  -»  0,  we  find 

C  =  i  fo r  all  time  (2.7) 

N 

Figure  2  demonstrates  that,  in  order  to  maintain  the  same  thinning  rate 

e  =  -  H/H  as  a  Newtonian  fluid  with  the  same  shear  viscosity  u  ,  one 
o 

would  have  to  impose  a  larger  pressure  P  on  this  viscoelastic  fluid. 

This  phenomenon  arises  not  from  viscoelasticity  per  se,  but  from  the 
elongational  viscosity  function  associated  with  this  particular  fluid. 

We  conclude  that  an  elongational  viscosity  of  the  type  often 
identified  with  polymer  solutions  exerts  a  retarding  influence  on  the 
thinning  rate  of  a  liquid  disk  subjected  to  a  uniform  pressure.  Thus, 
elongational  viscosity  may  play  an  important  role  in  control¬ 
ling  events  which  are  believed  to  lead  to  droplet  disruption. 

4 .  BUBBLE  COLLAP1E  AS  AN  ELONGATIONAL  FLOW 

13-15 

Earlier  studies  have  established  that  the  flow  field  in  a  liquid 

surrounding  a  spherical  collapsing  bubble  is  an  elongational  flow.  This 


16 


follows  directly  from  the  kinematics,  defined  by  the  continuity  equation 
in  spherical  coordinates.  Assuming  spherically  symmetric  collapse,  one 
integrates 


S_ 

dr 


Z 

r  v 

r 


=  0 


(28) 


to  find 


(29) 


where  R  -  dR/dt  is  the  rate  of  change  of  bubble  radius.  It  is  useful 
to  note,  at  this  point,  that  we  may  write 


A(t) 


-2  "  1  3Tr3 


(30) 
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The  components  of  y' 


11 


[o] 


are  given 


0 


0 


(33) 


Here,  r  and  r  '  are  the  radial  positions  of  a  material  element  of 
fluid  at  times  t  and  t',  respectively.  The  relationship  of  r  to 
r'  is  found  through  the  kinematics,  using 


3t 


A  (t  ) 


Integration  gives 


/  3  / 

r  (t  ) 


(34) 


(35) 


The  constant  of  integration  c'  is  found  using  the  fact  that  r'  =  r 
t '  -  t  .  It  follows  that 
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/3 

r 


3 

r 


+  (R 


r3» 


(36) 


and 


It  follows  further  that  the  components  of  y  ,  1  may  be  written  in  terms 

~l°J 

of  the  bubble  radius  R  upon  noting  that 


1  + 


R'3  -  R3 


.1/3 


(38) 


It  must  be  kept  in  mind  that  the  components  of  y  r  ,  are  functions  of 

~l°J 

r,  t'  ,  and  t. 


To  find  the  history  of  the  bubble  collapse,  R(t),  one  must  solve 
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the  Rayleigh  equation,  as  derived  in  our  previous  studies  in  the  form 


a 

p|rr  6  1  R2)  =  PR  -  P.  *  J 


[V 


T  1  dr 
—  r 


(39) 


R 


ZO 


As  in  the  previous  work1'5  we  may  reduce  thiB  to  the  form 

as 

,  /*♦  3  •  2\  f  T rr 

PG  -  Po  -  2y/R  -  p  RR  +  | R  J  =  -  2  I  - -  — —  dr  (40) 

R 

where  y  is  the  surface  tension,  and  -  P#  =  -  AP  is  the  pressure 

difference  between  the  bubble  interior  and  the  liquid  far  from  the  bubble. 

For  the  Newtonian  liquid,  for  which 

r  =  fi  A  i  1) 

we  find  R (t)  from  the  Rayleigh  equation  in  the  form 

-  AP  -  2y/R  -  p^RR  +  |  R2]  =  4p  R/R  (42) 

We  may  use  Eq.  42  as  an  aid  in  experimental  design,  and  we  will  return  to 
this  point  shortly. 

5.  EXPERIMENTAL  DESIGN 

The  design  concept  requires  that  a  small  gas  bubble  be  created  within 
the  liquid  of  interest.  An  excess  pressure  upon  the  liquid  will  then  cause  the 
bubble  to  collapse.  The  history  of  collapse,  R(t),  will  reflect  the  pxoperties 
o"  the  liquid,  and  will  be  described  by  a  solution  to  the  Rayleigh  equation, 

Eq.  40.  This  history  will  be  recorded  through  high  speed  photography. 

Prior  to  designing  an  experiment  it  is  necessary  to  establish  the 
working  range  of  parameters  that  will  be  used.  We  wish  to  use  the  observa¬ 


tion  of  bubble  collapse  to  infer  something  about  the  rheology  of  the  fluid 


surrounding  the  bubble. 

If  the  dynamics  is  too  rapid,  fluid  inertia  will  control  the  bubble 
response,  and  mask  the  role  of  rheology.  Thus,  it  is  essential  that  the  ex¬ 
periment  be  designed  so  as  to  minimize  ine  Tial  effects.  At  the  same  time, 
the  dynamics  must  be  fast  enough  so  that  th_^  high  deformation  rates  typical 
of  drop  shattering  are  simulated.  The  result  is  that  there  is  a  narrow 
"operating  window"  of  acceptable  parameters. 

We  can  use  the  Newtonian  form  of  the  Rayleigh  equation  to  quantify 
some  of  these  ideas. 

We  begin  with  the  assumption  that  the  collapse  is  so  fast  that  it  is 
adiabatic,  and  there  is  no  time  for  diffusion  of  gas  to  occur.  Hence  there 
is  a  constraint  between  pressure  and  volume.  Assuming  an  ideal  gas  we 
find 


T-,  ,rm 

P  V  =  constant 


(43) 


where  P  is  an  absolute  pressure,  and  where  m  is  the  heat  capacity 
Cj 

ratio 


m  =  C.  /C 
P  v 


(44) 


A  more  useful  form  is 


pcR 


3  m 


P  ,R3m 
o  C  j  o 


(45) 


where  P  is  the  initial  (absolute)  pressure  inside  the  bubble,  upon 

O  I 


dd 


initiation  of  collapse. 


Hence 


■£P  -  P  -  P  =  P  I  - 

*  G  oG  I  R 


P  -  P  +  P  I  —  -  1 

oG  “  oGil  R  I 


L  Po  +  PoG  \  R 


The  Rayleigh  equation  takes  the  form 


-iPo  +  PoG  T 


,  ••  3  •  2 

1  -  2y/R  -  P  RR  +  -  R 

2 


=  4  gR/R  (47) 


subject  to  initial  conditions 


R  =  R,  h  =  0  at  ts-  0 
o 


Note  that  for  t  <  0  ,  the  bubble  is  in  equilibrium  (R  =  R  )  at  a  pressure 
given  by  (setting  R  =  R  and  R  =  0  in  Eq.  47) 


£P  -  2y/R  =  0 
o  o 


P  +  2y/R  t  <  0 


Collapse  is  initiated  upon  incrementing  the  external  pressure  by 


an  amount  so  that 


AP  =  P  -  p  =  p  -  Zy)  R 
o  ®  oG  1  o 


(49) 


Upon  introducing  Eqs.  48  and  49  into  Eq.  4?  we  find 


-  P,  t  2y/R  +  IP  +  2y/R  ) 
1  o  l  06  o  / 


R 
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3  ft  2 

+  IR 


4p  R/R 


(50) 


We  may  first  inquire  into  the  total  extent  of  collapse.  For  long  time,  when 
R  approaches  zero,  we  find  from  Eq.  50  that 


P1  + 


2y/RQ  +  (p^  +  2y/RQj 


FR 

o 

R 


3m 


-1 


-  2y/R  =  0 
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For  the  sake  of  simple  calculation,  let  us  neglect  surface  tension 
in  Eq.  51  and  find 


R 

a 

T 


- 1  /  3m 


=  1  + 


P  +  P, 

*  l 


1  /3m 


(52) 


where  is  the  final  radius  of  the  collapsed  bubble.  We  see  that  the 

degree  of  collapse  depends  upon  the  ratio  of  absolute  pressures  before  and 
after  initiation  of  collapse. 
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If  we  take  a  reasonable  value  lor  ihe  heat  capacity  ratio  m  to  be 
m  =  1.4  for  diatomic  gases  we  can  prepare  a  simple  table  of  collapse 

ratios  as  a  function  of  pressure  ratios. 


Table.  Collapse  Ratios  as  a  Function  of 
Pressure  Ratios  (m  =1.4) 


According  to  this  estimate  (which  neglects  surface  tension)  a  pressure  ratio 
in  excess  of  ten  is  required  to  collapse  the  bubble  by  a  factor  of  two.  This 
is  a  sufficiently  accurate  estimate  for  the  purpose  of  design  of  the  experi¬ 
ment.  We  conclude  that  in  order  to  achieve  a  reasonable  degree  of  bubble 
collapse  (in  effect,  enough  data  to  be  meaningful)  we  will  have  to  use  pres¬ 
sure  ratios  of  order  ten. 

The  next  problem  to  deal  with  in  expe  rimental  design,  is  the  rate  of 
collapse,  since  this  will  determine  the  decree  of  need  for  high  speed  record 
ing  technology.  Again,  estimates  are  most  easily  made  using  Fq.  b'J  (i.e.  . 
assuming  Newtonian  behavior)  and  neglecting  surface  tension.  ,.\  simple, 
but  useful,  analytical  mode!  follows  if  we  nc.lcct  the  gas  pressure  inside  tin 


bubble.  This  provides  an  estimate  for  collapse  that  would  be  valid  for  short 
times,  and  small  extents  of  collapse.  Again,  we  note  that  the  goal  here  is 
the  development  of  simple  guides  to  experimental  design. 

If  we  go  back  to  the  Rayleigh  equation  in  the  form  of  Eq.  42,  and 
assume  that  AP  is  constant,  neglect  surface  tension,  and  neglect  inertial 
effects,  we  find 


A  useful  time  scale,  from  this  model,  is  the  time  to  collapse  to  1 /e  =  0.37 
of  the  initial  size.  (Note  that  under  constant  pressure,  ignoring  the  internal 
gas,  the  bubble  collapses  totally.  ) 

This  time  scale  is 

t  =  4p/AP  (54) 

For  a  10  poise  liquid  collapsing  under  a  100  psi  pressure  difference,  we  find 
t  =  6  x  10  ^s.  This  time  scale  is  too  short  for  existing  technology.  We 
note  that  liquids  of  interest  could  well  be  less  viscous  than  10  poise  (further 
shortening  the  time  scale).  This  result  suggests  that  it  will  be  necessary  to 
collapse  under  small  pressure  differences.  However,  the  previous  result 
on  extent  of  collapse  indicates  the  need  for  large  pressure  ratios.  It  may 
be  necessary  to  operate  at  sub-atmospheric  pressure,  so  that  1  +  P^  / Pm 
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remains  small.  For  example,  if  the  system  is  initially  at  0.  5  psia  (which 

requires  a  cold  sample  since  the  vapor  pressure  of  water,  for  example,  is 

not  much  less  than  this  at  room  temperature)  and  is  subjected  to  a  5  psia 

pressure,  so  that  AP  =  4.5  psia  while  P, /P  =  10,  then  we  might 

-4 

expect  a  measurable  bubble  collapse  with  a  time  scale  of  10  s.  It  is  clear 

that  we  will  be  dealing  with  very  high  speed  events. 

Another  parameter  that  can  be  controlled  is  the  bubble  size  R  . 

o 

If  is  too  large,  inertial  effects  dominate  the  collapse.  If  this  is  the 

case,  then  rheological  effects  are  unimportant,  and  it  is  not  possible  to  ex¬ 
tract  rheological  information  from  an  experiment  dominated  by  inertial  forces. 

Figure  3  shows  solutions  of  the  full  Rayleigh  equation  for  a  Newtonian 
liquid  for  bubbles  of  initial  radii  of  50  and  100  micrometers.  For  the  larger 
bubble,  and  the  lower  viscosity,  some  oscillation  of  the  bubble  occurs  before 
it  reaches  its  equilibrium  final  size.  This  is  an  inertial  effect.  For  10  and 
100  poise  liquids,  the  collapse  is  monotoruc.  Shown  on  Figure  3  is  the  nurely 
viscous  solution  (Eq.  53),  and  we  see  that  even  the  10  poise  liquid  deviates 
from  that  model  at  small  time,  again  demonstrating  some  degree  of  inertial 
resistance  to  collapse.  However,  the  100  poise  liquid  shows  purely  viscous 
collapse  over  most  of  its  history. 

By  contrast,  for  a  bubble  of  50  micrometers  radius,  no  oscillation 
occurs  in  the  5  poise  liquid,  and  nearly  viscous  behavioi  is  seen  at  10  poise. 

It  would  appear  from  these  simulations  that  it  will  be  necessary  to  use  bubbles 
of  50  micrometers  radius  or  less  in  order  to  minimize  inertial  effects. 
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6. 


THE  EXPERIMENTAL  SYSTEM 


6 .  1  Apparatus 

A  schematic  of  the  experimental  system  is  shown  in  figure  4,  and 
a  photograph  of  the  actual  apparatus  is  given  in  figure  5.  The  entire  system 
was  set  up  in  a  room  which  could  be  made  completely  dark,  for  photographic 
purposes.  The  experimental  system  may  be  conveniently  separated  into  four 
distinct  portions,  which  may  be  analyzed  individually.  The  numbers  in  paren¬ 
theses  refer  to  the  schematic  in  figure  4. 

6.  i.  i  Cell; 

The  heart  of  the  system  is  the  glass/copper  cell  (1)  with  the  platinum 
electrode  mounted  in  it.  Figures  6  and  7  are,  respectively,  a 
schematic  and  photograph  of  the  cell.  The  glass  is  standard 
glass  tubing  with  O.  D.  =  0.  25"  ,  I.  D.  =  0.020"  and 
L  =  1.5''.  The  copper  is  standard  copper  tubing  with  O.  D.  =  0.25", 
I.  D.  =  0.  185"  and  L  =  1.5".  The  glass  and  copper  tubes  are  ce¬ 
mented  together  using  Devcon1'®'  5-minute  epoxy.  The  acrylic 
supports  are  added  to  increase  strength  and  provide  a  way  to  mount 
the  cell  onto  the  microscope  stage. 

Figure  8  i6  a  schematic  of  the  platinum  electrode.  The  platinum 
wire  (O.  D.  -  100  pm)  is  fused  into  the  glass  sleeve  by  melting 
the  glass  around  the  wire.  The  top  ox  the  electrode  is  hhen  finely 
polished.  The  platinum  electrode  is  mounted  into  the  glass  cell 
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Glass  Cell 


using  a  high  temperature  epoxy.  Two  leads  come  off  the  cell 
to  the  pulse  generator  (12),  one  from  the  platinum  wire  at  the 
bottom  of  the  cell  and  the  other  from  the  copper  tubing  at  the 
top.  When  the  glass/copper  cell  is  filled  with  a  conducting 

liquid  (i.  e.  ,  the  test  solution)  the  electric  circuit  is  completed. 

(R) 

Bubbles  are  generated  using  a  Tektronix^  160A  Power  Supply, 

161  Pulse  Generator  and  162  Waveform  Generator  in  series  (12). 
Bubbles  of  various  sizes  can  be  made  by  varying  pulse  width 
and  magnitude.  The  electric  pulse  hydrolyzes  the  water  and 
forms  a  hydrogen  gas  bubble  which  then  rises  in  the  liquid. 

Optical  System: 

The  optical  system  is  a  self-contained  unit  consisting  of  a  Leits^ 

;Si 

microscope  (2,  18)  with  a  Wildv>  MPS-11  Shutterpiece  mechanism 
(15)  attached.  The  microscope  assembly  is  mounted  in  a  hori¬ 
zontal  position  (as  opposed  to  the  usual  vertical,  stand-up 
position)  on  a  wooden  platform  which  can  be  raised  and  lowered. 
This  allows  for  easy  alignment  with  the  rotary  drum  camera  (3). 
The  cell  is  fixed  to  the  microscope  stage,  which  moves  left  to 
right  and  forward/backward  to  allow  centering  and  focusing  of 
the  rising  gas  bubble.  The  microscope  optics  were  used  as 
delivered,  but  the  lighting  system  was  modified  for  photographic 
purposes.  The  12  V,  50  W  Halogen-Bellaphot*^  (OSRAM)  lamp 
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The  Wilti^  MPS- 1 1  Shutterpiece  (15)  was  used  to  control  ex¬ 
posure  times.  This  uses  an  electronic  shutter  mechanism  so 
this  signal  was  also  used  to  trigger  the  solenoid  valve  in  the 
pressure  system.  Variable  shutter  speeds  are  available  but 
1/15  sec  proved  to  be  a  convenient  speed  at  which  to  conduct 
the  experiments.  A  black  cardboard  baffle  (16)  was  made  to 
connect  the  shutter  mechanism  to  the  lens  of  the  rotary  drum 
camera.  This  prevented  extraneous  light  from  entering  into 
the  lens.  The  lens  of  the  drum  camera  was  completely  covered 
except  for  a  vertical  slit  approximately  1/16"  wide.  This 
collimated  light  source  was  essential  to  obtaining  good  contrast 
on  the  streak  photographs.  Kodak  247  5  Recording  film  (1000 
ASA)  was  found  to  provide  the  best  balance  of  high  speed  and  good 
contrast. 

6.1.4  Pressure  System: 

To  collapse  the  bubble,  which  is  at  atmospheric  pressure  in  the 
cell,  it  was  necessary  to  pressurize  the  system.  The  source  of 
the  pressure  was  a  compressed  air  cylinder  (11)  controlled  by 
a  Purox ^  CGA- 5 80  pressure  regulator  (10).  To  instantaneously 

pressurize  the  system  a  fast  acting  valve  was  needed.  A  Schroder 

( R) 

Bellows^- cyclone  three-way  direct  acting  solenoid  valve  (9)  w."s 
chosen.  (Schroder  Part  ff  74323-01  15,  Orifice  1/16/1/16,  t20V/ 
60  Hz,  brass  body.  )  The  solenoid  valve  was  activated  by  an 
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electronic  trigger  from  the  MPS-11  Shutterpiece  via  an  electronic 
relay  (GE  #SC-43/530T)  (8).  In  this  way  synchronization  of  the 


pressure  pulse  and  shutter  opening  was  made  possible.  There 
was  a  short  time  delay  between  shutter  opening  and  pressurization 
of  the  system  so  it  was  certain  that  the  event  would  not  occur  be¬ 


fore  the  shutter  had  opened.  Pressure  was  measured  using  a 

(R) 

Valid yne^  DPI  S  variable  reluctance  differential  pressure  trans¬ 


ducer  (5)  in  conjunction  with  a  Validyne 


(R) 


CD15  sine  wave  carrier 


demodulator  (6),  connected  to  a  digital  voltmeter  (Hewlett-Packard) 


(7).  Various  pressure  plates  (diaphragms)  are  available  for  the 
transducer,  so  that  a  wide  pressure  range  can  be  covered.  The 

pressure  hose  was  connected  directly  to  the  top  of  the  cell  (copper 

(g\ 

section)  using  a  hose  clamp  and  some  RTV^Silicone  rubber 
sealant.  The  cell  was  able  to  withstand  pressures  up  to  160  psig. 


•  ^  Experimental  Procedure. 

The  fluid  to  be  used  is  loaded  into  the  cell  using  a  disposable  pipet. 
Care  must  be  taken  to  make  sure  that  the  fluid  level  is  beyond  the  glass/copper 
interface,  because  contact  with  the  copper  is  necessary  to  complete  the  circuit 
for  the  bubble  hydrolysis.  The  pressure  tubing  is  then  fastened  to  the  top  of 
the  cell,  using  both  a  hose  clamp  and  RTV  Silicone  Rubber  sealant  (GE).  The 
entire  assembly  is  mounted  onto  the  moveable  stage  of  the  microscope.  The 
leads  from  the  pulse  generator  are  attached  to  the  cell;  one  to  the  platinum 
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wire  which  extends  out  the  bottom  of  the  cell  and  the  other  to  the  copper 
tubing  on  top. 

Next,  the  room  is  made  completely  dark  and  the  film  is  loaded  into 
the  drum  camera.  Each  experimental  run  uses  one  35  mm  roll  (36  exposure) 
of  the  Kodak  2475  recording  film.  The  drum  camera  is  equipped  with  a  three 
position  lever  (shut,  focus,  and  open),  so  that  after  the  film  is  loaded  the 
lever  is  moved  to  the  shut  position  and  the  lights  in  the  experimental  system 
can  be  put  on.  The  drum  camera  is  set  in  motion  and  its  speed  monitored 
using  the  Tektronix'^  storage  oscilloscope.  A  shutter  speed  of  1/15  sec  is 
chosen  for  the  experiment  so  the  drum  speed  must  be  kept  below  15  rev/sec. 
Typically,  experiments  are  run  at  fi  to  12  rev/sec  (125  to  80  m  sec/rev  ). 
The  significance  of  this  lime  with  respect  to  the  streak  photographs  will  be 
explained  later. 

Next,  the  pulse  generator  and  pressure  system  are  turned  on. 
Various  size  bubbles  in  the  range  of  50  to  500  (Am  can  be  made,  depending 
on  pulse  width  and  magnitude.  In  general,  for  the  purpose  of  photographic 
simplicity,  bubbles  in  the  200  to  300  pm  range  are  used.  Once  made,  the 
bubble  rises  slowly  up  the  cell  until  it  conies  into  viewing  range.  (Note: 

In  the  fluids  that  were  studied  in  these  experiments,  it  was  found  that  the 
upward  movement  of  the  bubble  is  much  "slower"  than  the  time  scale  of  the 
experiment,  so  for  all  practical  purposes  the  bubble  is  treated  as  stationary. 


For  less  viscous  fluids,  such  as  pure  water  for  instance,  this  can  be  a  big 
problem,  )  At  this  point  the  lever  on  the  camera  is  moved  to  the  focus 
position,  and  the  bubble  can  be  aligned  and  focused  while  one  looks  through 
the  viewing  port  at  the  back  of  the  drum  camera.  The  moveable  microscope 
stage  makes  It  easy  to  center  and  focus  the  drop.  With  everything  in  position 
the  shutter  is  closed  and  set  to  1/15  sec,  the  lever  on  the  camera  is  put  in 
tne  open  position  and  the  pressure  on  the  air  tank  regulator  is  set  to  the 
desired  level  (typically  60  to  160  psig).  Now,  the  shutter  is  opened,  which 
(riggers  the  solenoid  valve  through  the  electronic  relay,  and  opens  the  cell 
to  die  desired  pressure.  The  bubble  collapses  and  the  event  is  recorded  on 
film  in  the  form  of  a  "streak  photograph"  (which  will  be  discussed  later). 

The  rotation  speed  of  the  drum  cameri  is  obtained  from  the  storage  scope 
and  the  pressure  is  read  from  a  calibrated  digital  voltmeter  attached  to  the 
sine  wave  demodulator. 

Once  again,  the  room  is  made  completely  dark  and  the  film  is 
removed  from  the  drum  camera,  loaded  onto  a  35  mm  developing  reel  and 
placed  in  a  standard  light-tight  developing  tank..  The  film  is  "push-processed" 
by  developing  for  11  to  12  min.  in  Kodak  DK-50  developer  (recommended  de¬ 
veloping  time  is  5  to  6  min.  ).  The  end  result  is  a  streak  photograph  of  the 
bubble  collapse,  as  is  shown  in  figure  9. 

6 .  3  Solutions . 

The  solutions  that  were  used  in  this  study  were  all  water  soluble 
polymer  systems.  The  two  polymers  used  were  carboxymethylcellulose  (CMC) 
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(Hercules,  Inc.  CMC  7H3S,  experimental  sample)  and  polyacrylamide 
(PAA)  (Stein,  Hall  and  Co.  ,  Polyhall  2  95).  The  four  solutions  that  were 
used  in  this  study  are  CMC  1%;  PAA  0.25,  1.0  and  1.  5%  where  X%  is 
defined  as  wt.  polymer/wt.  H^O. 

The  rheological  properties  of  these  samples  were  obtained  using 

(R) 

a  Rheometrics^ Fluids  Rheometer  in  the  cone  and  plate  geometry  and  the 

* 

dynamic  oscillatory  mode.  Figure  10  shows  the  complex  viscosity  (17  ) 
for  the  four  samples.  In  the  limit  of  us  -*  0 ,  77  -*  T)q  ,  the  steady  state 

zero-shear  viscosity,  and  over  the  entire  frequency  range  the  behavior  of  T) 
parallels  that  of  T).  From  the  plots  of  complex  viscosity  versus  frequency 
several  characteristics  of  the  fluids  are  apparent.  First,  in  every  case  the 
PAA  solutions  are  more  highly  non-newtonian  than  the  CMC  solution.  Second, 
below  a)  =  10  rad/sec,  the  CMC  solution  behaves  essentially  as  a  Newtonian 
fluid.  Hence,  we  hope  to  be  able  to  use  this  CMC  solution  as  a  psuedo- 
newtonian  fluid  for  experimental  purposes. 

The  dynamic  storage  (elastic)  and  loss  (viscous)  moduli,  c'  and 
C"  respectively,  for  the  four  solutions  are  given  in  figurj  11.  Once  again, 
there  is  a  marked  difference  in  the  behavior  of  the  CMC  and  PAA  solutions. 

In  the  case  of  CMC,  the  viscous  forces  (G  )  are  always  substantially 
greater  than  the  elastic  forces  (G  ).  For  the  PAA  solutions  at  the  low 
frequencies,  the  viscous  forces  (G")  are  larger  than  the  elastic  forces 
(G/),  but  as  the  frequency  gets  higher  the  two  forces  become  equal  in 
magnitude  and  finally  the  elastic  forces  begin  to  dominate.  This  markedly 
different  behavior  between  the  PAA  and  CMC  might  allow  one  to  speculate 
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that  the  elastic  forces  would  be  a  factor  in  one  case  and  not  in  the  other. 

It  is  also  evident  that  the  frequency  (strain  rate)  range  of  the  experiment 
could  be  very  important  in  determining  fluid  behavior  because  of  the  strong 
frequency  dependent  nature  of  the  polymer  solutions. 

7.  CALIBRATION  AND  ANALYSIS  OF  EXPERIMENTAL  DATA 

7.  I  Magnification  Calibration. 

To  determine  the  total  magnification  in  the  optical  system  a  simple 
test  was  conducted.  A  thin  chromel  wire  (O.  D.  =  0.  0093")  was  inserted  into 
the  cell  when  it  was  filled  with  a  test  fluid.  A  series  of  still  photographs  were 
taken  with  the  rotary  drum  camera.  From  the  wire  dimensions  measured  on 
film  and  the  known  (measured)  diameter  of  the  wire,  a  magnification  factor 
was  determined.  The  magnification  was  found  to  be  55X. 

7 .  2  Determination  of  Collapse  Times. 

Since  there  is  no  readily  available  time  measure  (such  as  framing 
rate)  in  the  rotary  drum  camera,  a  method  to  calculate  the  time  scale  of  the 
event  had  to  be  devised.  What  was  done  was  to  measure  the  perimeter  of  the 
drum  where  the  film  is  placed  during  the  experiment.  This  gives  a  measure 
of  length  per  revolution.  From  the  storage  scope,  the  rotational  speed  is 

obtained.  Thus,  one  is  able  to  calculate  length  per  second.  Typical  time 

-3  -4 

scales  for  collapse  were  approximately  10  to  10  sec.  From  the 
streak  photographs  it  is  possible  to  calculate  the  "lag  time"  between  the 
shutter  opening  and  the  pressure  pulse  hitting  the  fluid.  This  was  found  to 
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be  approximate!/  15  msecs. 


8.  EXPERlh.  :NTAL  DATA 

Figure  12  is  a  photograph  of  a  single  bubble  taken  with  the  rotary 
drum  camera  when  it  is  not  rotating  and  the  slit  has  been  removed  from  the 
camera  lens.  This  bubble  is  approximately  200  pm  in  diameter.  The  white 
line  in  the  center  of  the  bubble  is  caused  by  the  diffraction  of  light  as  it  passeB 
through  the  bubble. 

If  one  starts  with  the  bubble  as  pictured  in  figure  12  and  covers 
the  lens  with  a  very  narrow  vertical  slit,  all  that  can  be  seen  is  a  slice 
through  the  bubble.  Therefore,  with  the  camera  rotating,  the  collapse  appears 
as  a  streak  of  continually  decreasing  width.  This  is  shown  by  the  "streak 
photograph"  given  in  figure  9.  Referring  to  that  photograph,  one  sees  that 
there  is  an  i.  itial  portion,  which  is  indicative  of  the  time  between  the  shutter 
opening  and  the  pressure  pulse  hitting  the  system.  Once  the  pressure  hits, 
the  collapse  begins  (small  oscillations  can  be  seen  initially)  and  the  streak 
continues  to  decrease  in  width  until  the  equilibrium  bubble  size  is  reached. 

The  pressure  remains  constant  throughout  the  run. 

Analysis  of  the  data  consists  of  measuring  the  initial  streak  width 
and  the  width  as  a  function  of  distance  until  no  more  change  is  observed. 

The  streak  width  is  correlated  to  the  bubble  diameter  (using  the  magnification 
factor  previously  determined)  and  the  distance  along  the  film  related  to  the 
time  of  collapse  (using  the  rotation  speed  and  length  of  film  per  rotation). 
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Figure  13  shows  data  from  two  bubble  collapse  runs  —  one  for 
the  relatively  inelastic  CMC  solution  and  one  for  the  more  elastic  PAA 
solution.  The  most  striking  observation  is  the  similarity  of  collapse  of 
the  two  fluids.  Not  obvious  from  the  figure  is  the  fact  that  the  collapse 
is  much  slower  than  would  be  expected  for  these  fluids.  Finally  we  note 
that  the  collapse  proceeds  to  a  level  below  that  predicted  by  Eq.  52,  based 
on  the  adiabatic  assumption.  Similar  results  are  observed  in  all  cases 
studied  to  date.  Additional  data  are  in  the  appendix. 

These  results  provide  clues  to  problems  in  the  experimental 
system  that  must  be  r.orrected. 

9.  EVALUATION  AND  DISCUSSION  OF  DATA 

Let  us  begin  with  the  positive  accomplishments,  and  then  examine 
the  problems  that  remain. 

We  have  designed,  constructed,  and  tested  a  system  for  high  speed 
photography  of  bubble  collapse.  While  improvements  in  the  photographic 
system  are  possible,  it  is  clear  that  accurate  and  reliable  data  are  possible. 
The  prototype  test  cell  is  not  made  of  glass  of  any  optical  quality.  Additional 
experience  with  lighting  and  focusing  techniques  will  lead  to  some  improve¬ 
ment  in  the  contrast  and  sharpness  of  the  streak  image.  It  Bhould  be 
possible  to  reduce  the  bubble  size  by  factors  of  two  to  three  and  still  obtain 
reliable  data, 
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The  technique  for  production  of  bubbles  in  aqueous  solutions 
works  well.  In  these  preliminary  tests  relatively  large  bubbles  have 
been  used,  of  about  100  pm  radius.  It  is  not  difficult  to  reduce  the 
bubble  radius  to  30  to  50  pm,  and  we  will  note  subsequently  that  this 
may  be  necessary. 

The  major  problem  lies  in  the  imposition  of  a  step-change  in 
pressure  upon  the  test  cell.  Superimposed  upon  figure  13  is  a  trace 
of  the  pressure  measured  in  the  test  cell.  It  would  appear  that  upon 
opening  the  solenoid  valve  there  is  a  relatively  rapid  "spike"  of 
pressure  followed  by  a  slow  rise.  Thus,  the  bubble  is  not  collapsing 
under  an  instantaneous  step-change  to  the  new  pressure.  To  date,  it 
has  not  been  possible  to  determine  if  the  problem  lies  in  the  solenoid 
valve  itself,  or  in  some  other  feature  of  the  design  of  the  test  cell/valve/ 
pressure  reservoir  system.  This  would  appear  to  be  the  major  problem 
to  be  dealt  with  at  this  point  in  the  development  of  the  experimental 
system.  Until  this  problem  is  solved,  data  obtained  to  date  cannot  be 
compared  to  any  theoretical  model  based  upon  solution  of  the  appropriate 
form  of  the  Rayleigh  equation. 
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10. 


RECOMME  ND  ATIONS 


10.  1  There  are  several  areas  where  improvements  could  be  made 

in  the  experimental  system  to  extend  its  usefulness  and  accuracy. 

10.  1.  1  Cell. 

A  new  cell  could  be  made  of  higher  quality  glass.  Special  opti¬ 
cal  glass  cells,  as  are  used  in  spectrophotometers,  are  available  com¬ 
mercially.  It  might  be  useful  to  make  the  cell  rectangular  in  order  to 
prevent  any  optical  distortion  that  may  occur  when  photographing  through 
a  circular  cell. 

10.1.2  Optical  System. 

It  appears  that  most  nf  the  photographic  problems  which  were 
encountered  earlier  in  the  project  have  been  resolved.  We  would  like  to 
obtain  sharper  photographs,  but  a  new  cell  should  help  considerably  in 
this  area.  The  keys  to  obtaining  high  contrast  photographs  are  the  very- 
high  intensity  light  and  the  use  of  the  slit  on  the  camera  lens  to  exclude 
stray  light. 

10.  1.  3  Pressure  System. 

This  is  the  area  where  the  biggest  improvements  could  be 
made.  The  solenoid  valve  used  in  our  experiments  is,  to  our  knowledge, 
the  fastest  opening  valve  available  for  this  pressure  range.  Unless  some 


improvement  in  this  area  of  technolgy  exists,  of  which  we  are  unaware, 
it  would  appear  that  a  solenoid  valve  w'iil  not  be  suitable  for  use  with 
low  visdosity  fluids  whose  collapse  times  are  of  the  order  examined  in 
our  studies.  The  current  technology  is  suited  to  "slower"  fluids,  i.  e.  , 
more  viscous  fluids,  or  to  collapse  under  smaller  driving  pressures 
(see  10.  2.  3  below). 

An  alternative  to  a  solenoid  valve  should  be  explored.  One 
possibility  would  be  to  separate  the  test  cell  and  the  pressure  reservoir 
by  a  thin  membrane  or  septum  that  could  be  punctured  with  a  mechanical 
device.  This  alternative  was  not  examined,  and  at  this  point  we  arc  un¬ 
able  to  estimate  the  character  of  the  pressure  pulse  that  such  a  system 
would  produce. 

We  conclude  that  the  limitations  of  current  technology  for 
creating  a  change  in  pressure  make  it  impossible  to  impose  a  suitable 
approximation  to  a  "step  change"  in  pressure  on  the  bubble.  However, 
the  pressure  pulse  can  be  measured,  and  the  theoretical  analysis  for 
bubble  collapse  could  be  suitably  modified  to  account  for  the  observed 
pressure  pulse. 

10.  2  In  addition  to  these  hardware  improvements,  several  changes 

in  procedure  would  lead  to  more  meaningful  results. 

1  0  -  2 .  !  Bubble  Si  v.v  . 


Large  bubbles  have  been  used  in  these  preliminary  tests  in 


order  to  minimize  photographic  problems.  A  redaction  by  a  factor  of 
two  to  three,  down  to  the  range  of  35  to  50  pm  radius,  is  possible  with 
the  current  design,  and  will  ensure  the  elimination  of  inertial  effects. 

10.  2.  2  Bubbles  in  Non- Aqueous  Solutions. 

At  the  present  time  the  bubble  production  system  can  only 
be  used  for  aqueous  materials.  To  extend  the  use  of  the  system  to 
organic  materials,  another  means  of  bubble  formation,  not  based  on 
hydrolysis,  is  necessary. 

A  new  technique  for  making  small  bubbles,  suitable  for  both 
aqueous  and  organic  solvents,  has  been  discovered.  Some  of  the  solu¬ 
tion  of  interest  is  placed  in  a  high  speed  blender  and  beaten  into  a  froth. 

A  sample  of  this  aerated  liquid  is  removed.  Most  of  the  bubbles,  because 
of  their  size,  will  rise  and  leave  the  liquid.  High  speed  agitation  of  a 
viscous  liquid  will  create  some  bubbles  oi  extremely  small  size.  The 
desired  range  is  50  to  100  microns  diameter.  These  small  bubbles  will 
remain  in  the  sample.  With  a  hypodermic,  a  very  small  volume  of  the 
sample  is  removed.  The  volume  is  such  as  to  contain  only  a  few  small 
bubbles.  This  sample  is  then  introduced  into  the  bottom  of  the  test  cell 
which  contains  the  quiescent  test  liquid.  When  one  of  the  small  bubbles 
rises  into  the  main  body  of  the  test  liquid,  the  collapse  experiment  is 
initiated.  Specific  work  remains  to  establish  an  optimum  procedure 
with  respect  to  agitation,  sample  sH  ,  etc.  A  modified  test  cell  will 
have  to  be  designed. 
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10.  2.  3 


Reduced  Pressures. 


Ultimately  one  will  wish  to  study  more  dilute,  lower  viscosity, 
solutions.  In  order  to  slow  down  the  collapse,  it  will  be  necessary  to  use 
smaller  driving  pressures.  However,  it  is  still  necessary  to  use  large 
pressure  ratios  in  order  to  achieve  significant  changes  in  bubble  size. 

This  will  necessitate  the  use  of  sub-atmospheric  pressure  prior  to  col¬ 
lapse,  followed  by  imposition  of  a  small  increment  of  pressure.  For 
example,  absolute  pressures  before  and  during  collapse  couid  be  o.  5  psia 
and  5  psia.  In  this  way  the  required  large  pressure  ratio  of  order  ten  is 
maintained,  while  the  pressure  difference  is  only  4.  5  psi.  One  possible 
advantage  of  this  is  the  likelihood  that  it  will  be  easier  to  produce  an 
"instantaneous"  step-change  in  pressure  at  lower  pressure  differences. 

Of  course,  the  lower  limit  on  pressure  is  the  vapor  pres¬ 
sure  of  the  solution.  However,  one  could  work  with  a  cold  solution, 

and  extrapolate  resulto  to  (higher)  temperatures  of  interest  through 

,  il 

the  use  of  classical  time-temperature  superposition  principles. 


54 


LITERATURE  CITED 


1.  Adelberg,  M.  Mean  Drop  Size  Resulting  from  the  Injection  of 
a  Liquid  Jet  into  a  High-Speed  Gas  Stream.  ALAA  J.  6.,  1143 
(1968). 

2.  Wilcox,  J.  D.  ,  R.  K.  June,  H.  A.  Brown,  Jr.  ,  and  R.  C. 

Kelly,  Jr.  The  Retardation  of  Drop  Breakup  in  High  Velocity 
Airstreams  by  Polymeric  Modifiers.  J.  Appl,  Polymer.  Sci.  JS, 
1  (1961). 

3.  Wang,  K-H.  ,  and  C.  Tien.  Atomization  and  Drop  Size  of 
Polymer  Solution,  IEC  Proc.  Des.  Dev.  _H,  169  (1972). 

4.  Gerber,  B.  V.,  and  A.  K.  Stuempfle.  TSD  004.  7- 11.  AD- 
AC25999/4GA.  A  New  Experimental  Technique  for  Studying 
the  Explosive  Comminution  of  Liquids.  Available  from  NTIS 
U.  S.  Dept.  Commerce.  1976  Army  Science  Conference 
Proceedings.  Vol.  I.  1976  UNCLASSIFIED  Report. 

5.  Simpkins,  P.  G.  ,  and  E.  L.  Bales.  Water-Drop  Response 
to  Sudden  Accelerations.  J.  Fluid  Mech,  55..  629  (1972). 

6.  Fishburn,  B.  D.  Boundary  Layer  Stripping  of  Liquid  Drops 
Fragmented  by  Taylor  Instability.  Acta  Astron.  2>  1267  (1974). 

7.  Hanson,  A.  R.  ,  E.  G.  Domich,  and  H.  S.  Adams.  Shock- 
Tube  Investigation  of  the  Breakup  of  Drops  by  Air  Blasts.  Phys. 
Fluids.  6,  1070  (1963). 

8.  Haas,  F.  C.  Stability  of  Droplets  Suddenly  Exposed  to  a  High 
Velocity  Gas  Stream.  A.  I.  Ch.  E.  Journal.  J_0,  920  (1964). 

9.  Harper,  E.  Y.  ,  G.  W.  Grube,  and  I-D.  Chang.  On  the 
Breakup  of  Accelerating  Liquid  Drops.  J.  Fluid  Mech.  52,  - 
565  (1972). 

10.  Ranger,  A,  A.  ,  and  J.  A,  Nicholls,  Aerodynamic  Shattering 
of  Liquid  Drops.  ALAA  Journal.  2-  285  (1969). 

11.  Bird,  R.  B.  ,  R.  C.  Armstrong,  and  O.  Hassager.  Dynamics 
of  Polymeric  Liquids.  Vols.  1  and  2.  Wiley,  NY.  1977. 


55 


12. 


Shir  jdkar,  P.  ,  and  S.  Middleman.  Lubrication  Flows  in 
Viscoelastic  Liquids:  1.  Squeezing  Flow  between  Approaching 
Parallel  Rigid  Planes.  J.  Rheology. 26^  xx  (1982). 

Pearson,  G.  H.  and  S.  Middleman.  Elongational  Flow  Behavior 
of  Viscoelastic  Liquids:  Part  I.  Modeling  of  Bubble  Collapse. 
AIChE  Journal.  23,  714  (1977). 

Pearson,  G,  H,  and  S.  Middleman.  Elongational  Flow  Behavior 
of  Viscoelastic  Liquids:  Part  II.  Definition  and  Measurement  of 
Apparent  Elongational  Viscosity.  AIChE  Journal.  II,  111  (1977) 

P®^rson,  G.  H.  and  S.  Middleman.  Elongational  Flow  Behavior 
of  Viscoelastic  Liquids:  Modeling  Bubble  Dynamics  with  Visco¬ 
elastic  Constitutive  Equations.  Rheol.  Acta.  17,  500  (1978). 


APPENDIX 


(tMcl)  ■*U,  u,e*lil*3  &  <me*  ;//» 


GLOSSARY 


2  , 

Bo  =  gPrQ/y 

g 

H 

M 


rr. 


P 


r 

o 

R 

r 

T 

t 

V 

v 

z 

a 


=  t/A 

=  (vr>  Vg,  vz) 

-  i  (-H) 

2  H 


Bond  number 
Acceleration  of  gravity 
Disk  thickness 
Memory  function 
Heat  capacity  ratio 

Increment  in  pressure  causing  collapse 

Pressure  inside  the  bubble 

External  pressure  on  the  fluid 

Pressure  distribution  v/ithin  fluid 

Drop  radius 

Radius  of  bubble 

Radial  coordinate 

Time  variable 
Volume  of  bubble 
Velocity  vector 
Axial  coordinate 


y 


Surface  tension 


Strain  tensor 


y 

A 

€  =  -  H/H 

o 

V 

X 

A  =  X  f 
o 

M- 

n 

p 

r 


Shear  rate 

Rate  C'f  deformation  tensor 

Nonnewtonian  viscosity 
Relaxation  Lime 

Newtonian  viscosity 
Dimensionless  pressure 
Fluid  density 
Stress  tensor 


66 


DISTRIBUTION  LIST  4  j 

Names 

Cop  1  as 

1 

Names  Copies  -I 

CHEMICAL  systems  laboratory 

-  A 

Commander  1 

ATTN: 

DRDAR-CLB 

1 

USASED,  Uc  *.  1  NSCON  . 

ATTN : 

DRDAR-CL3-C 

1 

ATTN:  1 AFM-SED-  Ml  1  .  J 

ATTN: 

DRDAR-CLB-PO 

1 

Fort  Mea  'e,  MO  20755  1 

ATTN: 

DRDAR-CL8-PS 

! 

j 

ATTN: 

DRDAR-CLB-R 

I 

DEPARTMENT  OF  THE  ARMY  i 

ATTN  : 

DRDAR-CLB-R  ( A ) 

1 

ATTN: 

DRQAR-CLB-R ( M ) 

1 

hqda  J 

ATTN: 

DRDAR-CLB-R(S) 

1 

ATTN:  DAMO-NCC  1  ! 

ATTN: 

QRDAR-CLB-T 

t 

ATTN:  DAMO-NC/COI.  Robinson  1 

ATTN: 

DRDAR-CIB-T  £ 

I 

WASH  DC  20310 

ATTN: 

DROAR-CLC-3 

! 

ATTN: 

DRDAR-CLC-C 

1 

Federal  Emergency  Management  Agency  - 

ATTN: 

DRDAR-CLC-E 

1 

Office  of  Research/NPP 

ATTN: 

DRDAR-CLF 

1 

ATTN:  Da  v  1  d  W  •  Be  risen  1 

ATTN: 

DRDAR-Ct J-P 

2 

Washington,  DC  20472 

ATTN  : 

DRDAR-CL J-L 

2 

ATTN: 

DROAR-CLJ-M 

! 

H?  CA  | 

ATTN: 

DROAR-CL  N 

1 

Office  of  the  Deputy  Chief  of  Staff  for 

ATTN: 

ORDAR-CLN-TE 

1 

Research,  Development  A  Acoulsltlon 

ATTN: 

DROAR-CL  T 

1 

ATTN;  DAMA-CSS-C  1 

ATTN: 

DRDAR-CLW-C 

1 

Washington,  DC  20310 

ATTN: 

DRDAR-CLW-P 

1 

ATTN: 

ORDAR-CLW-E 

1 

HQ  S 1 xth  US  Army 

ATTN: 

DRCAR-CLY-A 

1 

ATTN:  AFKC-OP-NSC  1 

ATTN- 

DRDAR-C  L Y -R 

4 

presidio  of  San  Francisco,  CA  94  129 

CCP  1  E S 

FOR  AUThOR(S) 

Co  mma  nder 

Research  Division  (CPC) 

6 

DARC0M,  STITEUR  j 

RECORD 

COPY:  DRCAR-CLS-A 

1 

ATTN;  DRXST-STI  1 

Box  48,  AD0  No*  York  09710  4 

DEPARTMENT  Of  DEFENSE 

Commander  - 

Defense  Technical  Information 

Center 

USASTCFEO 

ATTN: 

DT  IC-ODA-2 

2 

ATTN:  MAJ  Ml  keworth  1  a 

Cameron  Station,  Building  5 

APD  San  Franc'sco  96328  _ 

Alexandria,  VA  77  3  14 

i 

, 

HQ,  5th  Infantry  Division  (Mech) 

0 i rector 

ATTN:  Division  Chemical  Officer  1 

Defense  Intelligence  Agency 

Fort  Po  Ik,  LA  7  1459 

ATTN; 

D5-4G1 

1 

wash  1  n 

gto.i,  DC  20  30  1 

I 

67 

4 

i 

Comma  rider 

US  Army  Nuclear  A  Chemical  Agency 
ATTN:  MONA-WE  I 

7500  Back  I  I ck  fid.  Bldg  2073 
Springfield,  VA  22150 

Army  Research  Office 

ATTN:  DRXRO-Cfi  (Dr.  R.  Gh Ira-  del  II)  1 

P.0,  Box  12211 

Research  Triangle  Park,  NC  27709 
OFFICE  OF  THE  SURGEON  GENERAL 

Commander 

US  Army  Meclcal  Bioengineering  Research 
and  Development  Laboratory 
ATTN:  SGRD-UBD-AL,  Bldg  568  1 

Fort  Detrick,  Frederick,  MO  21701 

Commander 

USA  Medical  Research  Institute  of 
Chemical  Defense 

ATTN:  SGRD-UV-L  1 

*t>9'-doen  Proving  Ground,  MO  2101C 

US  ARMY  MATERIEL  DEVELOPMENT  AND 
READINESS  COMMAND 

Commander 

US  Army  Materiel  Development 
Readiness  Command 
ATTN:  DRCLOC 

ATTN;  DRCSF-P 
5001  Elsnnhowor  Ave 
Alexandria,  VA  22333 

PM  Smoke/Obscurants 
ATTN:  DRCPM-SFK-S  3 

Aberdeen  Proving  Ground,  MD  21005 


Dl rector 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-MP  1 

ATTN:  DRXSY-CA  (Mr.  Metz)  1 

Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Army  Missile  Command 
Redstone  Scientific  Information  Center 
ATTN:  DRSMI-RPR  (Documents)  1 

Redstone  Arsenal,  AL  35809 

Dl rector 

DARC0M  Field  Safety  Activity 

ATTN:  DRX0S-C  I 

Charlestown,  IN  47111 


Commander 

US  Army  Natick  Research  and  Development 
La  borator I es 


ATTN: 

DRDNA-0 

1 

ATTN; 

DRDNA-  1C 

1 

ATTN: 

DRDNA- ICAA 

1 

ATTN: 

DRDNA-  IV. 

1 

ATTN: 

DRDNA- ITF  (Dr.  Roy  W.  Roth) 

2 

Natick, 

MA  01760 

US  ARMY 

ARMAMENT  RESEARCH  AND 

DEVELOPMENT  COMMAND 

Commander 

US  Army 

Armament  Research  and 

Development  Command 

ATTN: 

DRDAR-  LCA-L 

1 

ATTN: 

DRDAR-LCU-CE 

1 

ATTN: 

DRDAR- NC  (COL  Lymn) 

3 

ATTN: 

DRDAR-SCM 

1 

ATTN: 

DRDAR-SCS 

1 

ATTN: 

DRDAR-TDC  (Dr.  D.  Gyorog) 

1 

ATTN: 

DRDAR- T  S  S 

2 

ATTN  : 

ORC 

1 

Do  vor. 

NJ 

Domnan  der 

US  Army  Foreign  Sc  I '  ■  co  A  Technology  Center 
ATTN:  DRX5T-MT3  1 

220  Seventh  St , ,  NL 
Char  I  o 1 1  e s v  I  I  I  o,  VA  2290  1 


and 

I 

1 


60 


US  Army  Armament  Research  and 
Development  Command 
Resident  Operations  Office 
ATTN:  DRDAR-TSE-OA  (Robert  Thresher) 
National  Space  Technology  laboratories 
NSTi  Station,  MS  39529 


Comma  ndant 
US  Ar  my  Missile  I 
and  School 
ATTN:  ATSK-CM 

ATTN:  ATSK-TME 

Redstone  Arsenal, 


Hun  1 1 Ions  Cen ' er 


AL  358C9 


Commander 

ARRAOCOM 

ATTN:  DRDAR-QAC-E 

Aberdeen  Proving  Ground,  MD  21010 


Commander 

US  Army  Logistics  Center 
1  ATTN:  ATCL-MG 

Fort  Lee,  VA  23801 


Comma  nder 

USA  Technical  Detachment 
US  Naval  EOD  Technology  Center 
Indian  Head,  MD  20640 

US  ARMY  ARMAMENT  MATERIEL  READINESS 
COMMAND 


Commanoant 

1  US  Army  Chemical  School 

ATTN:  ATZN-CM-C 

ATTN:  ATZN-CM-A0 

ATTN:  ATZN-CM-TPC 

Fort  McClel Ian,  AL  36205 


Commander 

US  Army  Armament  Materiel  Readiness  Command 


ATTN:  DRSAR-ASN 

ATTN:  DRSAR-IRW 
Rock  Island,  IL 


1 

1 

61299 


Comma  nder 
USA  ARRCOM 

ATTN:  CARTE  I 

Aberdeen  Proving  Ground,  MO  21010 


Comma  nder 

US  Army  Djgway  Proving  Ground 

ATTN:  Technical  Library  (Docu  Sect)  I 

Dugway,  UT  84022 


US  ARMY  TRAINING  4  DOCTRINE  COMMAND 


Commandant 

US  Army  Infantry  School 

ATTN:  CTD0,  CSD,  NBC  Branch  I 

fort  Bennlng,  GA  31905 


Commander 

USAAVNC 

ATTN:  ATZO-D-MS 

Fort  Rucker,  AL  36362 

Commander 

US  Army  Infantry  Center 
ATTN:  ATSM-CD-MS-C 

Fort  Bennlng,  GA  31905 

Commander 

US  Army  Infantry  Center 
Directorate  of  Plans  4  Training 
ATTN:  ATZB-DPT-PO-NBC 

Fort  Bennlng,  GA  31905 

Commander 

USA  Training  and  Doctrine  Command 

ATTN:  ATCD-.N 

Fort  Monroe,  VA  23651 


Commander 

US  Army  Ar  mor  Center 
ATTN:  ATZK-CO-  f-'s 

ATTN:  AT7K-PPT-P0-C 

Fort  Knox,  KY  40121 


Comma  nder 

USA  Combined  Arms  Canter  and 
Fort  Leavenworth 
ATTN:  ATZL-CAM-IM 
Fort  Leavenworth,  KS  66027 

US  ARMY  TEST  4  EVALUATION  COMMAND 

Commander 

US  Army  Test  4  Evaluation  Command 
ATTN:  ORSTE-CM-F 

ATTN:  DRSTE-CT-T 

Aberdeen  Proving  Ground,  MD  21005 
DEPARTMENT  OF  THE  NAVY 
Project  Manager 

Theatre  Nuclear  Warfare  Project  Office 
ATTN:  TN-09C 

Navy  Department 
Washington,  DC  20360 

Chief  of  Naval  Research 
ATTN :  Code  44  I 

800  N,  Quincy  Streot 
Arlington,  VA  22217 

Comma  nder 

Naval  Explosive  Ordnance  Disposal 
Technology  Center 
ATTN:  AC- 3 

Indian  Hoad,  Mn  20640 

Of  f  I  cer-  I  n-Ch  arge 
Marine  Corps  Detachment 
Naval  Explosive  Ordnance  Disposal 
Technology  Center 
Indian  Head,  MO  20640 

Comma  nder 

Naval  Surface  Weapons  Center 
Code  G5I 

L’ahlg-nn,  VA  22448 

^  n I e  f ,  Bureau  of  Mod. cine  4  Surgory 
Department  of  the  Navy 

ATTN:  MED  3033 

Washington,  DC  20372 


Commander 

Naval  Air  Development  Center 

ATTN:  Code  2012  (Dr,  Robert  Hetmbold)  I 

Warminster,  PA  18974 

US  MARINE  CORPS 

Commandant 

HQ,  US  Mar  I ne  Corps 

ATTN:  Code  LMW-50  1 

Washington,  DC  20380 

Commanding  General 
Marine  Corps  Development  and 
Education  Command 

ATTN:  Fire  Power  Division,  D09 1  1 

Quantlco,  VA  22134 

DEPARTMENT  OF  THE  AIR  FORCE 

AS0/AESD  , 

Wr  I  gh t-Pa t ter  son  A  F  B ,  OH  4  5433 

HQ  AFSC/SD2  , 

ATTN:  CPT  D.  Rledlgor 

Andrews  AFB ,  MC  20334 

HQ,  A  F  SC/S  CNE  ) 

Andrews  AFB,  M0  20334 

HQ,  A  F  SC/SGS  , 

Andrews  AFB,  00  20334 

HQ,  NORAD 

ATTN:  J-3TU  1 

Peterson  AFB,  CO  80914 

AFAMRL/HE 

ATTN;  Or,  Clyde  Replogqle  f 

Wr I ght-Patterson  AFB,  OH  45433 

HQ  AFTEC/TEL  , 

K  I  rt I  and  AFB,  NM  87117 

USAF  TAWC/THL  | 

Eg  I  I n  AFB,  FI  3254  2 


70 


AFATL/DLV 

Eglln  AF6,  FL  32542 


1 


OUTSIDE  AGENCI ES 


USAF  SC 

ATTN;  AD/YQ  1 

ATTN:  AD/YQO  (MAJ  Owens)  I 

Eglln  AFB,  FL  32542 

USAFSAM/VN 

Oeputy  for  Chemical  Defense 

ATTN:  Dr.  F.  Wesley  Baumgardner  I 

Brooks  AFB,  TX  78235 

AF AMR  L/TS 

ATTN;  COl  Johnson  I 

Wrlght-Patterson  AFB ,  OH  45433 


AMJ/RDTK 

ATTN:  LTC  T.  Kl ngery  1 

8rooks  AFB,  TX  78235 

AfO/RDSM  I 

Brooks  AFB ,  TX  78235 

AMD/RDSX  I 

Brooks  AFB,  TX  78235 

AO/XRO  I 

Eglln  AFB,  FL  32542 


Battel le,  Columbus  Laboratories 
ATTN:  TACTEC  1 

505  King  Avenue 
Co  I umbus ,  OH  4  3201 

Toxicology  Information  Center,  JH  652 
National  Research  Council  1 

2101  Constitution  Ave.  ,  NW 
Washington,  DC  20418 

US  Public  Health  Service 
Center  for  Disease  Control 

ATTN:  Lew  I  s  Webb,  Jr.  1 

Building  4,  Room  232 
Atlanta,  GA  30333 

D I  rector 

Central  Intelligence  Agency 

ATTN:  AMR/0RD/DD/S4T  1 

Washington,  DC  20505 

ADDITIONAL  ADDRESSEES 

Commandant 

Academy  of  Health  Sciences,  US  Army 
ATTN:  HSHA-CDH  1 

ATTN:  HSHA-IPM  2 

Fort  Sam  Houston,  TX  78234 


